Pion Interferometry for a Granular Source of Quark-Gluon Plasma Droplets 
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We examine the two-pion interferometry for a granular source of quark-gluon plasma droplets. 
The evolution of the droplets is described by relativistic hydrodynamics with an equation of state 
suggested by lattice gauge results. Pions are assumed to be emitted thermally from the droplets at 
the freeze-out configuration characterized by a freeze-out temperature Tj. We find that the HBT 
radius -R ou t decreases if the initial size of the droplets decreases. On the other hand, i? s id G depends 
on the droplet spatial distribution and is relatively independent of the droplet size. It increases 
with an increase in the width of the spatial distribution and the collective-expansion velocity of the 
droplets. As a result, the value of Rout can lie close to 7? s idc for a granular quark-gluon plasma 
source. The granular model of the emitting source may provide an explanation to the RHIC HBT 
puzzle and may lead to a new insight into the dynamics of the quark-gluon plasma phase transition. 

PACS numbers: 25.75.-q, 25.75.Nq, 25.75.Gz 
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Recent experimental pion HBT measurements at 
RHIC give the ratio of R out /i? s idc ~ 1 QLyL which is con- 
trary to many theoretical expectations "0, 0, EL IE IE 0@ • 
This RHIC HBT puzzle hints that the pion emitting time 
may be very short [9|, [lOj • Various models have been put 
forth to explain the HBT puzzle HJ El El El EE E3 

In this paper, we propose a simple model of a g ranular 
source of quark-gluon plasma droplets [TH fTst Ha. l2fll2lT | 
to explain the puzzle. The possible occurrence of a granu- 
lar structure of droplets during a first-order QCD phase 
transition was discussed by Witten and examined 
by many authors |23l 1 2 iL 123 . 1 2(1 l27j . We assume that 
pions are emitted thermally from hydrodynamically ex- 
panding droplets at a freeze-out temperature Tt and use 
relativistic hydrodynamics to describe the evolution of 
quark-gluon plasma droplets with an equation of state 
suggested by QCD lattice gauge results H |M El HI 
The two-pion correlation function can then be calculated 
after knowing the hydrodynamical solution pfl I2H l3l| . 
As the average freeze-out time is approximately propor- 
tional to the initial droplet size, we would like to see 
whether a granular source with many smaller droplets 
and their corresponding smaller freeze-out times will lead 
to .R ou t close to i? s ido- 

We study a quark-gluon plasma with no net-baryon 
content and use a equation of state of the fluid in terms 
of the entropy density function s(T) by H El El EE] 
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Here dQ and dn are respectively the degrees of freedom 
in the quark-gluon plasma phase and the hadronic phase, 
T c « 160 MeV is the transition temperature, s c is the 
entropy density at T c , and AT (between and 0.1T C ) is 
the width of the transition [lEg- In this paper, we take 
d Q = 37, d H = 3, T c = 160 MeV as in Ref. HISl, and 
take AT = 0.05T C . 



We shall make the approximation that the hydrody- 
namical solution for many independent droplets can be 
obtained by superposing the hydrodynamical solution of 
a single droplet. It suffices to focus attention on the hy- 
drodynamics of a single droplet. Knowing the entropy 
density s(T), one can get the pressure p, energy density 
e, and the velocity of sound cq in the droplet with the 
following equations as in Ref. H E^] , 



cff's(T'), 



e = Ts-p, 



dp 
~de 



(2) 



The energy momentum tensor of a thermalized fluid 
cell in the center-of-mass frame of the droplet is [3T. I32I . I33] ! 



T» v {x) = [e(x) + p(x)] u"(x)u u (x) - p{x)g^ , (3) 

where x is the space-time coordinate, u M = 7(1,1?) is the 
4-velocity of the cell, and g^ v is the metric tensor. With 
the local conservation of energy and momentum, one can 
follow Rischke and Gyulassy and get the equations for 
spherical geometry as H I2E] 



d t E + d r [{E + p)v] = -F, 
d t M + d r {Mv+p) = -G, 



where E = T 00 , M = T 0r , 
2v 
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We assume the initial conditions as 



e (0,r) = «{ '°' r ; rd ' v(0,r) 
1 0, r > rd, v ' 



0, r < r d , 

1, r > r d , 



(4) 
(5) 
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where eo = 1.875T c s c H |2£| is the initial energy den- 
sity of the droplets, and r d is the initial droplet radius. 
Using the Harten-Lax-van Leer-Einfeldt (HLLE) scheme 
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expressed as [37 



FIG. 1: (a) Temperature profile and (b) velocity profile for 
the droplet at t n = 3nAr<2,A = 0.99. (c) Isotherms for the 
droplet. 



|3l I28ll34l l35| and the relation of p = p(e) obtained from 
Eqs. Q and J2J), one can get the solution of the hydrody- 
namical equations for F = G = 0. One then obtains the 
solution for Eqs. (0J and JSJ by using the Sod's operator 
splitting method [j, l28l |36| . The grid spacing for the 
HLLE scheme is taken to be Ax = O.Olr^, and the time 
step for the HLLE scheme and Sod's method corrector 
is At = 0.99Aa; HHU. Figure 1(a) and 1(b) show the 
temperature and velocity profiles of the droplet. Figure 
1(c) gives the isotherms for the droplet. 

The two-particle Bose-Einstein correlation function is 
defined as the ratio of the two-particle momentum dis- 
tribution P(fci,fc 2 ) to the product of the single-particle 
momentum distribution P{ki)P{k 2 ). For a chaotic pion- 
emitting source, P(h) (i — 1,2) and P(ki,k 2 ) can be 



P(k u k 2 ) = \^{k 1 ,k 2 ;X 1 ,X 2 



(8) 



(9) 



where A(ki,Xi) is the magnitude of the amplitude for 
emitting a pion with 4-momentum ki — (ki,Ei) in the 
laboratory frame at Xi and is given by the Bose-Einstein 
distribution in the local rest frame of the source point. 
$(fci, k 2 ; Xi,X 2 ) is the two-pion wave function. Neglect- 
ing the absorption of the emitted pions by other droplets, 
<f>(/ci,fc 2 ; Xi 1 X 2 ) is simply 



§{k x ,k 2 ;X x ,x 2 ) 
1 



^/2 



A(k 1 ,X 1 )A{k 2 ,X 2 )e 



iki -Xx+ik-z-Xz 



iki -X 2 +ik 2 X\ 



(10) 



A{k 1 ,X 2 )A(k 2 ,X 1 )e 



Using the components of "out" and "side" |3, [7|, |38|, |39| of 
the relative momentum of the two pions, q = |ki — ka|, 
as variables, we can construct the correlation function 
C(gout, g s ido) from P{ki,k 2 ) and P{ki)P{k 2 ) by summing 
over ki and k 2 for each ((?out, 9sidc) bin. The HBT radius 
i? out and i? s ido can then be extracted by fitting the calcu- 
lated correlation function C(g ut, Qside) with the following 
parametrized correlation function 



C(<fout, fedc) = 1 + Ae q "uti 



'I ryA 

-^side^sid 



(11) 



The explicit procedure for calculating the two-pion cor- 
relation function is as follows. 

Step 1: select the emission points of the two pions 
randomly on the space-time freeze-out surfaces of the 
droplets, and get their space-time coordinate X\ and X 2 
in the laboratory frame. 

Step 2: generate the momenta k[ and k' 2 of the two 
pions in local frame according the Bose-Einstein distri- 
bution characterized by the temperature T/, and obtain 
their momenta k\ and k 2 in the laboratory frame by 
Lorentz transforms. 

Step 3: calculate [E[ /E x ] [E' 2 /E 2 ] for P(k l )P(k 2 ) and 
[E'JE^E'z/E^cosKki-faKXi-Xi)] for P(k u k 2 ), and 
accumulate them in the corresponding (<? ut,9sidc) bin. 

Step 4: repeat steps 1 through 4 many times to get the 
correlation function within a certain accuracy. 

We first examine the two-pion correlation function for 
a singlet droplet source. By fitting the two-dimension 
correlation function C(<7 ut, <7sidc) obtained with the 
above steps with Eq. (|llfl . we get the parameters i? ut, 
i? s idc, and A simultaneously. In our calculations, the 
transverse momenta of the pions are integrated over. The 
average transverse momenta of the pions in our fitting 
samples are 307, 329, and 386 MeV for T f =0.8T C , 0.65T C , 
and 0.5T C , respectively. The reason for a larger average 
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FIG. 2: The two-pion interferometry results for the one- 
droplet source as a function of the initial droplet radius r^. 
The symbols o, •. and * are for the freeze-out temperatures 
T f = 0.80T C , Tf = 0.65T C , and T f = 0.50T C , respectively. 



transverse momentum to associate with a smaller freeze- 
out temperature is due to the larger average expansion 
velocity in the case of a smaller freeze-out temperature. 
Figure 2(a), 2(b), 2(c), and 2(d) show the HBT results 
-Rout, -Rsidc, -Rout /-Rside, and A as a function of the ini- 
tial radius of the droplet, for the freeze-out tempera- 
tures Tf = 0.80T C (symbol o), T f = 0.65T C (symbol •), 
and Tf = 0.50T C (symbol *). It can be seen that the 
HBT radii i? ut and i? s ide increase linearly with r^, but 
the ratio -R ut/^?sidc is about 3 within the errors. The 
radius -R s ide reflects the spatial size of the source and 
the radius i? ou t is related to the lifetime of the source 
HL HE HH ■ From the hydrodynamical solution in fig- 
ure 1(c), both the average freeze-out time and freeze-out 
radial distance increase with Td for different Tf. As a 
consequence, -R ou t /-Rsidc is insensitive to the values Td 
and Tf. The value of -R ou t/ -Rsidc ~ 3 for a single droplet 
is however much larger than the observed values Q, • 
In our calculations, we did not including resonances in 
the hadronic phase. If we take the hypothetical case of 
dQ/dn — 3 to include a resonance gas in the hadronic 
phase, as discussed by Rischke and Gyulassy 0, we find 
that the ratio of i? ou t /-Rside is about 2.75. It is still much 
larger than the observed values 0,0- In Fig. 1 (d), the 
values of A for = 4 fm are larger than unity. This is 
a non-Gaussian correlation function effect, and the effect 
is larger for a wider correlation function, corresponding 
to a smaller source. 

As the average freeze-out time is proportional to the 
initial radius of the droplet, the freeze-out time and -R ou t 
decreases if the initial radius of the droplet decreases. On 
the other hand, -R s ide increases if the width of the droplet 
spatial distribution increases. A variation of the droplet 
size and the width of droplet spatial distribution can re- 
sult in -Rout nearly equal to -R s idc ■ Accordingly, we calcu- 



FIG. 3: Two-pion HBT results for the granular source of the 
droplets. The symbols o, •, and * are for the initial droplet 
radius ra —2.0, 1.5, and 1.0 fm, respectively. 



late next the two-pion correlation function for a Gaussian 
distribution source of Nd droplets. The spatial center- 
of-mass coordinates X<j of the droplets are assumed to 
obey a static Gaussian distribution^ exp(— X^/2-Rq) 
[HiiiJ. Figure 3(a), 3(b), 3(c), and 3(d) give the 
HBT -Rout, -Rside, -Rout /-Rsidc, and A as a function of the 
number of droplets Nd for different values of r^. In this 
calculation, we take Tf — 0.65T C and Rg — 5.0 fm. The 
symbols o, •, and * correspond to = 2.0 fm, Td = 1.5 
fm, and Td = 1.0 fm, respectively. It can be seen that the 
radii -R ou t and -Rside have a slowly increasing tendency 
as Nd increases but their ratio -R ou t / -Rside is almost inde- 
pendent of Nd- -Rout decreases as decreases but i? s ide 
is relatively independent of r^. Consequently, the ra- 
tio -Rout /-Rside decreases when decreases. The ratio 
-Rout /-Rside for Td — 1.5 fm is about 1.15 which is much 
smaller than the result of about 3 for the single droplet 
source. 

Finally, to study the effect of additional collective ex- 
pansion of the droplets, we calculate the two-pion cor- 
relation function for an expanding source. The initial 
distribution of the droplets is the same Gaussian distri- 
bution as the above granular source, but the droplets 
are assumed to expand collectively with a constant ra- 
dial velocity Vd after the initial time, in addition to their 
hydrodynamical expansion. Figure 4(a), 4(b), 4(c), and 
4(d) give the HBT R out , R sidc , -R ou t /-Rside , and A for dif- 
ferent values of Vd- In this calculation, we take the initial 
radius of the droplets to be Td = 1.5 fm, the other param- 
eters are the same as the above calculations for a static 
granular source. The symbols o, •, and * correspond 
to Vd =0, 0.3, and 0.6, respectively. The results in Figs. 
4(a) and 4(b) show that -R S ide increases more rapidly with 
the droplet collective expansion velocity Vd than i? out . A 
radial expansion will increase the transverse size of the 
granular source and -R s idc- On the other hand, -R ou t mea- 
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FIG. 4: Two-pion HBT results for the granular source with 
a radial collective expansion of the droplets. The symbols o, 
• , and * are for the expansion velocity Vd =0, 0.3, and 0.6, 
respectively. 



sures the source life time and the spatial extension where 
the two pions are emitted with nearly parallel and equal 
momenta, and the additional radial boost modifies only 
slightly the spatial separation of these points for most 
cases. As a result, i? ut does not increase as rapidly as 
i? s ido and i?out/-Rsidc is smaller at large Vd than at zero Vd 
(see Fig 1(c)). The ratio -R s ide/-Rout is of order 1 which 
is close to the observed value 0, Q . 

In summary, we propose a simple model of granular 
source of quark-gluon plasma droplets to examine the 
HBT interferometry data. The droplets evolve hydro- 
dynamically and pions are emitted thermally from the 
droplets at the freeze-out configuration characterized by 
a freeze-out temperature Tf. As the average freeze-out 



time is proportional to the radius of the droplet, smaller 
droplet size allows pions to be emitted within a shorter 
time and the life-time of the source decreases, leading to 
a smaller HBT radius i? ut- On the other hand, the HBT 
radius i? s ide depends on the width of the spatial distribu- 
tion of the droplets and is insensitive to the initial size of 
the droplets. The ratio of R out to i? s ide decreases signif- 
icantly if the emitted source is granular in nature. Fur- 
thermore, i? s idc increases with the collective-expansion 
velocity of the droplets more rapidly than i? ou t. The 
ratio i?out/-R s ido is about 1.15 — 0.88 for the collective- 
expansion velocity of the droplets from zero to 0.6, for 
the granular source with a Gaussian initial radius 5 fm 
and a droplet initial radius 1.5 fm. This i?out/-Rside ratio 
is close to the experimental values 0, ■ The granular 
model of quark-gluon plasma may provide a possible ex- 
planation to the RHIC HBT puzzle. It may also lead to a 
new insight into the dynamics of the quark-gluon plasma 
phase transition as the formation of a granular structure 
is expected to occur in a first-order QCD phase transition 

mm mm mis 

In order to bring out the most important features, we 
have neglected the multiple scattering effects on HBT 
interferometry [3ll Eol l4l| , and have not considered how 
the granular nature of the plasma may arise from detailed 
phase transition dynamics |2^, [2(| [23 • The sizes of the 
droplets in a collision can also have a distribution. Future 
refinements of the present model to take into account 
these effects on R ou t/R s idc will be of great interest. 
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